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ABSTRACT 

We analyze the velocity residuals of 551 carbon stars relative to a rotating-disk model of the inner 
^ 70deg^ of the Large Magellanic Cloud (LMC). We find that the great majority of the stars in this 
sample are best fit as being due to two different populations, a young disk population containing 20% 
of the stars with a velocity dispersion of 8kms^^, and an old disk containing the remaining stars with 
a velocity dispersion of 22kms~"'^. The young disk population has a metallicity ~ 0.25 dex higher than 
the old disk. 

With less certainty, the data also suggest at the 2cr level that there may be a third kinematically distinct 
population that is moving towards us at 30 km/sec relative to the LMC, consistent with measurements 
of 21 cm velocities. If real, this population contains about 7% of the carbon stars in the sample. It could 
be a feature in the disk of the LMC or it could be tidal debris in the foreground or background. If it is 
tidal debris, this population could account for some or all of the microlensing events observed towards 
the LMC. 
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1. INTRODUCTION 

Carbon stars are an important tracer of the kinematics 
of the disk of the Large Magellanic Cloud (LMC) . Kunkel 
et al. (1997) have analyzed the velocities of carbon s tars in 
the outer LMC disk. Hardy, Schommer, & Suntzeff ( |l999|) 
have measured the radial velocities of 551 carbon stars in 
the inner ~ 70 deg^ of the LMC and fit these velocities to 
a disk model. Here we focus on the residuals to this disk 
solution in order to isolate different kinematic components 
of the LMC carbon star population. 

The Milky Way disk has a multi-components structure 
that may be describable as a kinematically cold thin disk 
and a hotter thick disk as originally advocated by Gilmore 
& Reid (1983), or may comprise a continuum of struc- 
tures of increasing thickness (Norris & Ryan 1991). The 
LMC provides a unique laboratory to study the kinematic 
substructure of a quite different galaxy that also has disk 
kinematics. By studying this substructure, we can eventu- 
ally learn about the relations between formation history, 
disk heating, and enrichment in a non-Milky Way setting. 

Schommer et al. (1992) and Hughes, Wood, & Reid 
(1991) have presented data suggested that older com- 
ponents of the LMC disk are kinematically hotter than 
younger components, although these studies are limited 
by poor statistics. In this paper, we analyze a much larger 
radial-velocity sample to extract a substantially more de- 
tailed picture of the LMC's disk substructure. From both 
conceptual and practical standpoints, the analysis of the 
Carbon star radial velocities i s best divided into two steps. 
In the first step Hardy et al. (1999) fit for the global prop- 



erties of the disk including its projected rotation curve and 
its transverse velocity. Here we apply the second step and 
examine the residuals to that fits in order to extract in- 
formation about the kinematic structure of the disk. This 
study yields unambiguous evidence that the LMC disk, like 
the Milky Way disk, has a multi-component structure. We 
go on to show that, just as with the Milky Way, the colder 
disk component is more metal rich than the hotter one. 

In addition to determining the structure of the LMC 
disk, we also search for a non-disk component. One of the 
motivations for this research is the microlensing conun- 
drum. At present ~ 20 microlensing events towar ds the 
Magellanic clo uds h ave been analyzed (Alcock et al. 1997; 
Lasserre et al. 2000). If these microlensing events are due 
to halo objects, or Machos, then the detected Machos 
make up 10 — 30% of the mass of the halo. All obvious 
astrophysical candidates for halo microlensing have severe 
problems (e.g. Graff, Freese, Walker & Pinsonneault 1999) 
An alternative hypothesis is that the microlensing event s 
are due to lenses within the LMC (Wu [1994 Sahu |1994| ). 
However, if these lenses are virialized, they must have a 
large velocity dispersion (Gould 1995). In that case, we 
should see this population in the carbon star velocities, 
unless the carbon stars do not trace the lens population 
(Aubourg et al. 1999). 

Another possibility is that the observed microlensing is 
due to an unvirialized foreground or backgro und p opula- 
tion of lenses, such as a tidal stre amer (Zhao 199^; Zarit 



sky & Lin (1997); Zaritsky et al. |l999|) . In thiscase, we 



would expect the velocities of the lenses to be different 
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from those of LMC stars. Again, we should see this popu- 
lation in the carbon star velocities, unless the carbon stars 
do not trace the lens population or unless, by coincidence, 
the lens population has the same radial velocity as the 
main LMC population. We find that the data provide evi- 
dence at the 2 a level for additional velocity structure that 
could be due to an unvirialized foreground or background 
population. While this detection cannot be regarded as 
compelling, the problem of explaining the observed mi- 
crolensing events by other routes has proven so difficult 
that this proposed solution should be given serious con- 
sideration: our marginal detection should be checked by a 
much larger radial-velocity study. 

2. THE DATA 



Hardy et al. ( 1999| ) obtained radial velocities v for 551 
carbon stars in 35 fields, each about 0.25 deg^ scattered 
more or less uniformly over the inner 70deg^ of the LMC. 
The m easur ement errors are typically ~ lkms~^. Hardy 
et al. ( 1999 ) fit these velocities to a planar, inclined disk 
with a circular velocity that is allowed to vary in 5 bins. 
Table 1 shows a summary of the parameters for th e so- 
lution used in t his pa per; see Hardy et al. ( 1999 ) and 
Schommer et al. ( 1992 ) for details and descriptions of the 
rotation curve parameters and other possible fits. The 
fit adopted here is basically a solid body rotation model 
(constant dV/dr) out to 3.5 degrees, a flat rotation curve 
beyond that (3.5-5.5°), a slightly twisting line of nodes 
(0 in Table 1), an overall dispersion around the fit [a) 
which is characteristic of an intermediate to old disk pop- 
ulation, and an orbital transverse motion consistent with 
the proper motion measures of the LMC (e.g., Kroupa & 
Bastian 1997). The solution simultaneously fits for the 
transverse velocity of the LMC vj^ since this gives rise to 
a gradient in radial velocities across the face of the LMC 
with respect to angular position, Vw — v^. In this paper 
we primarily use the residuals to this fit. Aw, (§ 3 and § 
4.1) but also make use of the heliocentric radial velocities, 

V (§ 4.2). 

Table 1. Rotation Curve Parameters 

^ sys dV/dr ^circ 

50 km/s 21.5 km/s/kpc 75 km/s 

< e(PA) > 



-20° 



G Vtr 

18-22 km/s 250 km/sec 



3. DETECTION OF TWO POPULATIONS 

A histogram of the residuals Av is shown in Figure |[ 
We attempt to represent these residuals as various sums 
of Gaussians of the form 



N,, 



P(A«) = ^-=^exp 



2n(Tj 



{Av - Av,) 
2a2 



(1) 



subject to the constraint J2i ^i = 551. Here n is the num- 
ber of Gaussian components, and for each component i, 
Ni is the number of stars, Avi is the mean residual veloc- 
ity, and ai is the dispersion. We fit the velocity residuals 
to these functional forms by adjusting the parameters to 
maximize the log likelihood estimator. 
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lnL = ^ln[P(A 



Vk)]. 



(2) 



This is equivalent to a x^ minimization measurement in 
the Poisson limit of infinitely small bin size. Probabili- 
ties can be inferred from the log likelihood estimator by 
comparing likelihoods to the solution with maximum like- 
lihood and using the relation 



Ax' 



-2Alni. 



(3) 



Figure H shows fits to the (unbinned) residuals using a 
single Gaussian (with two free parameters) and a double 
Gaussian. In the latter fit, we impose the physically plau- 
sible additional constraint Awi — Aw2, so there are a total 
of 4 free parameters. The double Gaussian solution has 
20% of the stars in a thin disk population with a velocity 
dispersion of 8kms~^, and the remaining 80% of the stars 
in a thicker disk population with a velocity dispersion of 
22kms^^. The improvement is Ax^ — 20 for the addition 
of two degrees of freedom, i.e. a statistical significance of 
1 - exp(-AxV2) ~ 1 - lO-''-^. Thus, LMC carbon stars 
are better represented as two populations than one. How- 
ever, this does not prove that we have detected two distinct 
populations. It could also be that there are a continuum 
of populations with a range of dispersions from below 8 
to above 22kms~^. Nevertheless, for clarity of discussion, 
we will refer to two discrete populations. 

3.1. Metallicity of the two populations 



Costa & Frogel ( p96| ) (CF) pubhshed RI photometry of 
888 LMC carbon stars and 204 with infrared (JHK) pho- 
tometry. Within this sample, 103 of the stars that have 
infr ared p hotometry had velocities measured by Hardy et 
al. ( 19991 ). CF showed that the infrared colors differ be- 
tween samples of carbon stars from the Milky Way, the 
LMC, and the SMC. The carbon stars in the three galax- 
ies can be fit by 



{J - H)o = 0.62{H - K)o + C 



(4) 



with C, ~ {0.72, 0.67, 0.60} respectiv ely for the Galaxy, the 
LMC, and the SMC. Cohen et al. ( |l98l| ) suggested that 
this shift in colors is due to a metallicity related blanketing 
effect, in which case (^ can be used as a metallicity indica- 
tor. As can be seen in Figure 5 of CF, there is substantial 
scatter in the color-color relations compared to the dif- 
ferences among the three galaxies. Thus, this metallicity 
indicator cannot reliably determine the metallicity of an 
individual carbon star: it should be used only as a statis- 
tical estimator for stellar populations. 

Even though the metallicities of carbon stars in the 
three galaxies are unknown, if we assume that [Fe/H] ~ 
{0, —0.4, —0.8} for the three galaxies, we can make a rough 
calibration of this metallicity indicator: 



5[Fe/H] «6.7(5C. 



(5) 



k=l 



This relation should be taken only as rough estimate. 
However, one can be more confident of the relative or- 
der of the metallicities of carbon stars in the three galax- 
ies, and hence C can robustly distinguish between a high- 
metallicity population and a low-metallicity population. 

We find that the metallicity indicator ( is different for 
high velocity-residual stars than for low velocity stars. 
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Specifically, for stars with \Av\ < 10 km/sec, we find 
C = 0.678 ± 0.007 while for \Av\ > 10 km/sec, we have 
C = 0.662 ± 0.005. These two values of C are different at 
the 93% confidence level. However, since most of the "low 
velocity" stars chosen this way are actually from the more 
numerous thick-disk velocity sample, dividing up the sam- 
ple in this way is not the best way to measure the metal- 
licity difference. To isolate the thin and thick disks, we 
modify equation (|l|) to read 



A^, 



^(^-)-E:7^-p 



27rcr,- 



{Av - Av,) 



21 



2a 



exp 



(C-C.) 



21 



2(72 



where Q is the mean value of C for each population and 
(T|^ = 0.044 is the observed dispersion of C in the sample for 
the 103 stars with velocities and infrared data. Note that 
for stars without infrared data, the last term is simply set 
to umty. We then find (i = 0.663±0.04, C2 = 0.700±0.16, 
and C2 — Ci — 0.037 ± 0.017, i.e. a 2cr difference, which 
corresponds to A[Fe/H] ~ 0.25. 

Given the combination of different velocities and differ- 
ent metallicities, we claim that we have detected either 
two different disks within the LMC representing different 
ages of stellar populations or a continuous distribution of 
disk populations with a range of ages. In either case, the 
younger populations have higher metallicity and lower ve- 
locity dispersion. 

3.2. No virialized lenses 



Gould (1995) showed that for microlensing within a viri- 
alized disk, the microlensing optical depth is 



in the LMC, - 30kms-i (Hughes, Wood & Reid 1991, 
Schommer et al. 1992). Thus, we can show that the bulk 
of disk heating occurred during the Carbon star formation 
epoch. 



SEARCH FOR A KINEMATICALLY DISTINCT 
POPULATION 



The analysis of Gould (1995) only applies to virialized 
populations. It is still possible that an unvirialized popu- 
lation of stars could be causing microlensing. Such a popu- 
lation might be a streamer of stellar material pulled out by 
tidal interactions between the LMC and t he Mil ky Way, 
or between the LMC and the SMC (Zhao |1998| ). Zarit- 
sky & Lin (1997) claimed that they may have seen such a 
streamer in LMC clump giants. This paper caused numer- 
ous counter-argu ments w hich are summarized and debated 
in (Zaritsky et al. 



199^) 
ul: 



Ibata, Lewis & Beaulieu (|1998) examined the velocities 



of 40 clump giants in the LMC of which 24 were candidate 
foregr ound stars according to the criteria of Zaritsky & Lin 
(199'/). Ibata et al. (1998) found no difference in the mean 
velocities of the candidate foreground stars and the other 
clump stars and concluded that these stars did not form a 
separate kinematic population from the LMC. Zaritsky et 



al. ( 1999| ) confirmed the results of Ibata et al. (1998) using 
a much larger sample of 190 candidate foreground clump 
stars. However, the carbon-star sample that we analyze 
here is potentially more sensitive to the presence of tidal 
streamers than either of these two clump-star samples, in 
part because it is larger (551 stars) and in part because 
the velocity errors are much smaller (~ lkms~^). 



T = 2 — TT- sec i 



(7) 



where i is the angle of inclination of the disk with respect 
to the line of sight, 30 - 40° in the case of the LMC. In 
the case of the carbon stars, the total velocity dispersion is 
21 kms~^ and thus the optical depth due a virialised stellar 
population traced by the carbon stars is < 2 x 10~*, much 
smaller than that measured by the MACHO experiment 
(Alcock et al. |l997| ) of l.2toi x lO""^. Thus, the virial- 
ized population traced by carbon stars cannot account for 
microlensing. However, a virialized population too old to 
be traced by carbon stars would not be seen in our data 
(Aubourg et al. 1999). 

3.3. Conclusion 

We have explicitly assumed that hotter, more metal 
poor population is older than the younger, metal rich 
population in analogy with the Milky Way, even though 
the LMC may have a different disk heating mechanism 
than the Milky Way. The age-velocity dispersion relation 
has b een confirmed previously by H ugues, Wood & Reid 
and Schommer et al. (1992). Since we detect a 



( |l99il ) 



metallicity difference based on our infrared colors within 
this population, we also determine that some noticeable 
metal enrichment occured during the Carbon star forma- 
tion epoch. 

The velocity dispersion of the thick disk component, 
22kms~-'^, is much higher than the thin disk, and is close 
to the velocity dispersion of the oldest objects measured 



4.1. Search for third population in disk- fit residuals 

We search the data for a non- virialized, kinematically 
distinct population (KDP) in two different ways. First, 
we fit the residuals to the disk solution to the sum of three 
Gaussians, two representing the LMC, and one for the 
KDP. That is, we apply equation (||) with n = 3. We find 
a solution which is somewhat better than the two Gaus- 
sian fit, Ax^ = 8 for a change of 4 degrees of freedom. The 
off-center KDP peak is found to be moving towards us at 
27kms"^ relative to the bulk of the LMC and to contain 
63 stars, about 10% of the total. Thus, the data suggest 
that there may be a KDP, but at a statistically weak level 
of confidence. A Monte Carlo simulation was performed 
to verify the statistical confidence (details of which are de- 
scribed in § 4.3) which showed that this third bump is only 
present at the 75% confidence level. The fit to the third 
bump is shown in Fig 0. 



4.2. Search for a third population in velocities 

In the model considered in the previous section, the 
KDP stars have a common motion relative to the LMC. 
Possibly, the KDP stars are moving steadily away from 
the LMC disk, or are not associated with the LMC disk. 
The KDP should be seen in the original heliocentric radial 
velocities v better than it is seen in the disk-fit residuals 
Av. We therefore fit the data to a functions of the form 
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PKA«)=Eti^exp 



{Av-Av,y 
2^ 



exp 



.(C-COl 



2a 



Nu 



V27r(T 



■ exp 



exp 



(v-{vKup+A^e^+Ayeyf) 



(C-Ckdp)'' 



(8) 



where {0x,9y) is its angular position on the sky, and A^ 
and Ay are planar coefhcients for the heliocentric velocity 
distribution of the KDP. This equation is similar to equa- 
tion (o) but we have replaced the Au in the KDP terms by 
w, i.e., we fit to the heliocentric rather than the residual 
velocities. The origen of our x-y coordinate system is at 
a = 5''21",(5 = -69°f7', with X increasing to the east 
and Y to the north. 

Initially, we set A^ — Ay — Q, so that there are same 
number of degrees of freedom as in the three-Gaussian fit 
to the residuals. We find no solution here that has a lower 
X^ than the two-Gaussian solution, implying that there is 
no evidence for the existence of a third population having 
a common heliocentric velocity outside the LMC disk. 

We therefore repeat the search, but allow A^ and Ay 
to vary as free parameters. We find that the likelihood is 
then maximized at very low values of the velocity disper- 
sion ctkdp ^ lkms~^. We reject these soultions as un- 
physical, and note that our fitting routines may have been 
falsely attracted to them as results of inevitable Poisson 
noise. 

We then find a solution with 39 stars in the KDP 
with vkdp = 16.4kms~^, A^ = 2.6kms^^ dcg^^. Ay ~ 
4.9kms~-'^deg~-^, ctkdp = 5kms^^, and Ckdp = 0.673. 
Relative to the two-Gaussian solution, this KDP solution 
has Ax^ = 16 for 6 additional parameters. 

Figure shows the residuals of the LMC stars with re- 
spect to the KDP. The KDP is shown as the strong peak 
of points around residual 0. Other small peaks are due to 
the clumped distribution of our stars in angle, and are not 
significant. 

There are not enough stars in the KDP to significantly 
determine if the KDP covers the entire face of the LMC 
or has a patchy distribution. 

4.3. Monte Carlo 

While the probability that any randomly chosen plane 
will come within ^ 5 kms~^ of a significant fraction of our 
sample stars is small (and well represented by the x^ test), 
there are a large number of independent planes that can 
be compared to the data. To obtain a more accurate as- 
sessment of the statistical significance of this detection, we 
perform a set of Monte Carlo simulations. In each simula- 
tion, we draw velocities randomly from the two-Gaussian 
distribution of disk residuals found in § 3. We then search 
for a KDP in the resulting heliocentric velocities in the 
same way we did for the actual data in § 4.1 and § 4.2. In 
order to make the simulations tractable, we ignore metal- 
licity information. This simplification is justified by the 
fact that the metallicity of the KDP measured in § 4.2 
is not significantly different from the "young disk" com- 
ponent. If metallicity is ignored then the external-plane 
solution shows an improvement of Ax^ = 14 for 5 addi- 
tional parameters, which is formally significant at the 98% 
level. However, we find that out of 407 simulations, there 



is Ax^ > 14 in 26 cases. Hence, our detection is significant 
only at the 94% level, roughly equivalent to 2cr. 

4.4. Evidence of the KDP in Other LMC components 

Given the intriguing signal we see in the C star veloci- 
ties, but also the marginal level of significance, it is worth 
exploring other possible signs of the KDP. One such tracer 
i s the 21cm gas emissio n, ma pped, e.g., by Luks & Rohlfs 
( |1992| ), and Kim et al. ( |199S| ). Luks & Rohlfs note that a 
lower velocity component ( "L-component" ) contains about 
19% of the HI gas in the LMC, is separated from the main 
velocity component by ^--^30 km/s. Although Kim et al. 
(199^) do not specifically comment on such a component 



in their paper based on higher spatial resolution HI imag- 
ing, a similar signal seems evident in their position- velocity 
maps (e.g.. Figs. 7a and 7b in their paper) at RA 05:37 - 
05:47 and DEC -30 to -120 arcmin. The standard inter- 
pretation of this substructure in gas is that it is due to 
hydrodynamic effects on gas within the LMC disk. How- 
ever, the correlation of the gas velocity "L-component" 
with the stellar KDP suggests that the gas may be outside 
the LMC disk. 

An intriguing but somewhat more ambiguous signature 
may be e vident in the CH star velocities of Cowley & 
Hartwick (1991). Velocities for a sample of ~80 CH stars 
show a low velocity asymmetric tail, consistent with a com- 
ponent at ~20 km/sec lower systematic velocity. Cowley 
& Hartwick (1991) even suggest that one explanation of 
this population is that it is a result of an earlier violent 
tidal encounter between the LMC-SMC system and the 
Milky Way. The small sample statistics and asymmetric 
spatial distribution of these stars make a more detailed 
exploration difficult. 

5. MICROLENSING INTERPRETATION 

We may have detected a kinematically distinct popula- 
tion of carbon stars in the direction of the LMC. If real, 
this population could be either a structure within the LMC 
disk or tidal debris that is well separated from the disk and 
hence either in front of or behind the LMC. If it is well 
separated from the LMC, then it would give rise to mi- 
crolensing: either it would be in front of the LMC and so 
would act as lenses, or it would be behind the LMC and 
would act as sources. 

The microlensing optical depth due to a thin sheet of 
stellar matter with density Si and the LMC with density 
S2 separated by a distance D which is small compared to 
the distance from the Sun to the LMC is: 



47rG^ E1E2 

TKDP = — t;-L>- 



Si 



(9) 



The distance between the two sheets, D, cannot be de- 
termined from velocity data alone. However, since the two 
sheets must have similar velocities, the tidal tail cannot 
be a random interloper in the halo, but must be some- 
how related to the LMC. Lacking further information, we 
make the somewhat ad hoc assumption that the material 
in the tidal tail has been moving away from the LMC at a 
constant velocity of 30kms~^ since close tidal encounter 
between th e LM C and the SMC, 200 Myr ago (Gardiner 
& Noguchi 1996). In that case, we have 



D ~ WKDP X 200 Myr ~ 5 kpc. 



(10) 
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In fact, it is likely that the foreground object has had its 
velocity substantially changed by gravitational interaction 
with the LMC, and to a lesser extent, the SMC and the 
Milky Way, so this calculation only indicates that the ob- 
ject could have moved several kpc from the LMC in the 
past 200 Myr. All the results of this section will hold if 
the object is several kpc either in front of or behind the 
LMC. 

The total surface mass density. Si + S2, can be esti- 
mated from the observed surface brightness of the LMC, 
which is i? '^ 21.2 mag arcsec"^ (De Vaucouleurs 1957) 
near the center. If we assume a mass to light ratio of 3 
(in solar units), this corresponds to a total surface mass 
density of 300 Mq pc'^. 

It is possible that the surface densities of the disk and 
KDP populations are not traced by carbon stars. Still, 
lacking further information, we estimate the optical depth 
by setting S1/S2 — 39/(551 — 39) according to the solution 
of S 4.2, we obtain 



is 100 kms ^, in which case the typical mass of a lens is 



6 X 10"^ 



D 



5 kpc' 



(11) 



This optical depth is substantially larger than the opti- 
cal depth due to a virialized disk population traced by 
the carbon stars (< 2 x 10~*). It is consistent with the 
value observed by the Macho collaboration (Alcock et al. 
2000). There could be more tidal material which we have 
not found in this search because its velocity is by chance 
too close to the velocity of the LMC, and which would 
raise the optical depth. If D were greater than 5 kpc, then 
tkdp would rise proportionately. 

The transverse motion of such a population with respect 
to the LMC is probably 70 kms~^, the circular orbital 
velocity of the LMC. To calculate the typical transverse 
velocity in a microlensing event, this velocity should be 
added in quadrature to all the other sources of transverse 
velocity. The stars in the LMC are orbiting about the 
LMC center with a tran sverse motion of 70 kms^^ 
kpc (Kunkel et al. 1997| ; Hardy et al 



1999) 



at 4 

The LMC 

system has a transverse velocit y with respect to the Sun 
of some 250 kms~^ (Hardy et al. 1999) which will translate 
to a projected transverse motion of 25 kms~^ (at 5 kpc 
from the LMC). Adding these velocities in quadrature, the 
derived typical transverse velocity of a microlensing event 



M-0.13Af, 







D 



5 kpc 



(12) 



This is significantly below the mean mass of stars in the 
neighborhood of the Sun (e.g. Gould, Bahcall, & Flynn 
1997), but the LMC may have a different mass function. 
However, it is important to recognize that if D is made 
larger so as to account for more of the optical depth, then 
the mean mass is driven lower 



M~ 0.075 Af, 







1 X 10- 



(13) 



6. CONCLUSION 



We report two primary new results, one with high sta- 
tistical confidence, one which is shakier, but perhaps more 
interesting if true. We show that Carbon stars in the LMC 
are divided into a hot and cold population, with a clear dif- 
ference in metallicity between the two populations. Thus, 
we show that the epoch of LMC disk heating had to occur 
during the Carbon Star formation epoch. 

We also show with less confidence the existence of a 
third population, outside the LMC. If this population is 
real, it suggests that some fraction of the Carbon stars 
in the LMC are not in the disk, and thus could explain 
microlensing events. Although at present the statistical 
significance of this detection is not enviable, this result is 
still the best extant solution to the microlensing conun- 
drum. The microlensing conundrum poses such a difficult 
problem that several extreme explanations have been pro- 
posed including mirror matter and cosmological popula- 
tions of population III white dwarfs (Graff, Freese, Walker 
& Pinnsoneault 1999, and references therin). The kinc- 
matically distinct population is unique amongst these ex- 
planations in that it in not only allowed by the data, but 
even supported by the data at the 95% confidence level, 
and requires no modifications to the standard models of 
Particle Physics or Cosmology. We thus present it as the 
strongest explanation of LMC microlensing. 

Work at Ohio State was supported in part by grant AST 
97-27520 from the NSF. 
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6 Velocity structure of LMC 

FIGURE CAPTIONS 

Fig. 1. — The residuals of the carbon_stars with respect to the LMC disk fit. The grey Une is the best fit single gaussian. The black line is 
the best fit two gaussian model of eq. (|l|) . 

Fig. 2. — A fit to the residuals with three gaussians. Although the third peak is not significant in the fit to the residual, it is shown to be 
statistically significant when searched for in velocities, and shows the location of the KDP. 



Fig. 3. — The residuals of the stellar velocities with respect to the KDP. The KDP stands out as a strong peak near residual = 0. 
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